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> ' Abstract 

t> : 

We review hitherto attempts to look at the multiparticle production processes from the In- 
formation Theory point of view (both in its extensive and nonextensive versions). 
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: 1 Why Information Theory? 

^ . The idea of using Information Theory (IT) to analyze the multiparticle production data is 
very old one and can be traced back to 1J. This work attempted to establish what the 
experimental data of that time are telling us. The problem was serious because there was 
a number of theoretical models based on apparently completely disparate physical assump- 
tions, all of which were claiming to provide fairly good description of data and therefore 
were in fierce competition between themselves. The working hypothesis of [I] was that ex- 
perimental data contain only limited amount of information, which was common to all of 
them. Models considered differed therefore in some other aspects, which, from the point of 
experimental data considered, were however irrelevant. As it turned out, this information 
was that: (i) not all available energy sfs is used for production of secondaries (i.e., existence 
of inelasticity) and (ii) that transverse momenta od produced particles are strongly damped 
(apparent one-dimensionality of the relevant phase space). After closer scrutiny of models 
it was revealed that, indeed, all of them have these features build in (some in explicit some 
in very implicit way). From the point of view of data all these models were simply equiv- 
alent and their differences in what concerns their other physical assumptions were simply 
nonexistent. 
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2 Basics of IT approach 



This result could be obtained only by quantifying the notion of information, i.e., by resorting 
to IT with all its mathematical machinery. Referring to [2] for list of the relevant references 
we shall briefly sketch its basics. To do this one has first to introduce some measure of 
information contained in a given probability distribution pi. It is usually provided Shannon 
information entropy defined as: 

S = -EiPikipi. (1) 

Actually the choice of entropy used as measure of information reflects our a priori knowledge 
about hadronizing system. Out of many possibilities jH] we shall present here results obtained 
using the so called Tsallis entropy [I], 

S q = -^-ZUl-pl). (2) 

It is characterized by the so called nonextensivity parameter q, which quantifies deviation 
from the usual additivity of the entropy, namely that for S q , for two independent systems A 
and B 

S q (A+B) = S g A + S q B + (1 — q)S q AS q B- (3) 

Notice that in the limit q — > 1 one recovers the previous form of Boltzmann-Gibbs-Shannon 
entropy ((TJ). 

The IT tells us that the available data should be described by using only a truly minimal 
amount of information in order to avoid any unfounded and unnecessary assumptions. This 
information is provided by a finite number of observables {Fk, k — 1, . . . , n} of some physical 
quantities obtained by means of Pi and defined as: 

n 

(F k ) = J2PiF£\ (4) 

i=i 

when Shannon entropy is used and as: 

n 

{F k ) q = Y.M q ^\ (5) 
i=i 

when one uses Tsallis entropy instead. With this information one is then looking for (nor- 
malized) probability distribution {p{\ which contains only information provided by {Fk} 
and nothing more, i.e., which contain minimal information. As minimal information means 
maximal corresponding information entropy one is therefore looking for {pi} maximizing 
this entropy subjected to the to constraints given by eqs.(@]) or (jSJ) which account for our a 
priori knowledge of the process under consideration. As a result one gets the most probable 
and least biased distribution describing these data, which is not influenced by anything else 
besides the available information. In case of Shannon information entropy it is 1 



Pi = — exp 



fc=i 



(6) 



1 Notice that using the entropic measure S = Y^i [Pi^ n Pi T (1 i Pi)hi(l ± Pi)] (which, however, has 
nothing to do with IT) would result instead in Bose-Einstein and Fermi-Dirac formulas: Pi — (1/Z) ■ 
[exp[/3(ei — fi)] =F 1] , where f3 and /j, are obtained solving two constraint equations given, respectively, by 
energy and number of particles conservation [5] . It must be also stressed that the final functional form of 
Pi depends also on the functional form of the constraint functions Fk(xi). For example, F(x) cx ln(x) and 
ln(l — x) type constrains lead to pi cx xf(l — Xi) 13 distributions. 
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(where Z is obtained from the normalization condition J27=iPi = 1) whereas in case of Tsallis 
entropy it is 

k=l 

,(«) 



M 1 



= PI = ex P<? 



(7) 

where Z q is obtained from the normalization condition Y%=\ pT' = 1 an d where 



cxp 



A 



l-(l-g) 



A 
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(8) 



3 Examples of applications 



In Fig. ^ we provide some selected examples of application of IT to describe single particle 
distributions in different reactions. In general fits are good and prefer values of the nonex- 
tensivity parameter q > 1 (except in the attempt to fit e + e~ data, cf. last panel of Fig. 
^ where q < 1 is preferred and even then one cannot reproduce all features of data) and 
therefore its meaning here is worth of comment. As has been shown in ^0] it is given by 
fluctuations in the parameter 1/A in the exponential distribution of the form ~ exp(— x/A), 
namely 2 : 



\ 2 \ / i \ 2 



q = q T = 1 ± J f— 2 • (9) 



When applied to pr distributions (like in the upper-right panel of Fig. HI cf. also [T2] and 
references therein) it can be therefore connected with fluctuations of what is usually regarded 
in thermodynamical descriptions of collisions as being the "temperature" of the hadronizing 
system. In what concerns rapidity distributions (the rest of Fig. [T]) as was argued in , 

q = ql = i + p (io) 

where k is parameter characterizing (together with mean multiplicity n) the so called Neg- 
ative Binomial multiplicity distribution P^sin). This is so because, as was shown in |13j . 

Pnbd{u) = J o dn— = T{1 + n)m ■ (7 + 1)fc+ „> (11) 

where 

k , 1 a 2 (n) a 2 (n) 1 . . 

i.e., Pnbd{ji) arises from Poisson distribution by fluctuating its mean multiplicity n using 
gamma distribution. In general qL dominates in the collision process (being of the order of 
qL ~ 1.2 in comparison to qx ~ 1.02). One also observes tendency that qr is smaller for bigger 
hadronization systems [HI E2! what agrees with suggestion jO] that reflecting fluctuations of 
temperature qr = 1 + 1/ C where C is the heat capacity of the system and as such is expected 
to grow with the collision volume. 

2 Although in 10 fluctuations were assumed to be described by gamma function this result is general and 
lead to introduction of idea of super statistics, cf., jll| . 
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4 Summary 



The question arises: what is the advantage of the IT method? To answer it let us first 
notice that in examples presented in Fig. ^IT was represented by most probable and least 
biased distribution ((HJ) describing allocation of given (known) number of secondaries in a 
(longitudinal) phase space defined by a given (or assumed, by using parameter of inelasticity 
K) available energy. It means then that IE] in (JBJ) one has only one term, k = 1, with 
Fi being energy of secondary under consideration and the lagrange multiplier Ai = ft, being 
the inverse of "temperature" (understood in the sense mentioned before). This formula 
is apparently identical with formula used by statistical models of hadronization, however, 
here both ft and Z are not free parameters (see JB| for discussion and references) but are 
obtained from the constraint equations (here energy conservation and normalization) 3 . The 
parameters fitted is the energy available for hadronization (i.e., inelasticity parameter K, 
cf., E], in some cases like e + e~ collisions and some AA collisions they are fixed by the 
requirement of experiment) and parameter q, which as we argue, defines the amount of 
dynamical, intrinsic fluctuations present in the hadronizing system. In case that data cannot 
be fitted by this method we should add some other constraints (as in pA case [8J, or turn 
to the true dynamical description (as is probably the case with e + e~ collisions). 

Let us end with remark that IT does not solve our dynamical problems. On the other 
hand it is the only approach which allows us to select a minimal number of indispensable 
hypothesis (assumptions) needed to reproduce experimental data under consideration. In 
this approach any new hypothesis are allowed only when discrepancy with some new (or 
additional) experimental results occur. The choice of the form of information entropy (here 
represented by parameter q) offers additional flexibility because, as was stressed here, q sum- 
marizes many possible dynamical effects (out of which we have stressed here fluctuations 4 ). 
Therefore assumptions tested by using methods of IT can serve as ideal starting point to 
build any dynamical model of hadronization process. 

Two of us (OU and GW) would like to acknowledge support obtained from the Bogolyubov- 
Infeld program in JINR. Partial support of the Polish State Committee for Scientific Research 
(KBN) (grant 621/E-78/SPUB/CERN/P-03/DZ4/99 (GW)) is also acknowledged. 

References 

[1] Y.-A.Chao, Nucl. Phys. B40 (1972) 475. 

[2] G.Wilk and Z.Wlodarczyk, Phys. Rev. D43 (1991) 794. 

[3] F. Topsce, Physica A340 (2004) 11. 

[4] C. Tsallis, in Nonextensive Statistical Mechanics and its Applications, S.Abe and 
Y.Okamoto (Eds.), Lecture Notes in Physics LPN560, Springer (2000), in Physica A340 
(2004) 1 and Physica A344 (2004) 718, and references therein. 

3 Similar situation is for the formula Q. Actually, in both cases normalization Z and Z q can be exchanged 
for the requirement of properly and exactly reproducing the multiplicity of secondaries produced in a given 
event. 

4 It can be argued that it summarizes also effects of correlations, especially those arising from production 
of resonances • Recently question of connecting q with other forms of correlations has been discussed 

in HS|. 



4 



A.M. Teweldeberhan, H.G. Miller and R. Tegen, Int. J. Mod. Phys. E12 (2003) 395. 

F.S. Navarra, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, Phys. Rev. D67 (2003) 114002. 

F.S. Navarra, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, Physica A340 (2004) 467. 

O.V. Utyuzh, GWilk and Z. Wlodarczyk, Multiparticle production processes from the 
Information Theory point of view, hep-ph/0503048, to be published in Acta Phys. Hung. 
A (HIP) (2005). 

F. S. Navarra, O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, Physica A344 (2004) 568. 

G. Wilk and Z.Wlodarczyk, Phys. Rev. Lett. 84 (2000) 2770; Chaos, Solitons and Frac- 
tals 13 (2002) 581 and Physica A305 (2002) 227. 

C. Beck and E.G.D. Cohen, Physica A322 (2003) 267. 

M. Biyajima, M. Kaneyama, T. Mizoguchi and G. Wilk, Eur. Phys. J. C40 (2005) 243. 

P. Carruthers and C.C. Shih, Int. J. Mod. Phys. A4 (1989) 5587. 

O.V. Utyuzh, G. Wilk and Z. Wlodarczyk, J. Phys. G26 (2000) L39. 

C.Tsallis, M. Gell-Mann and Y.Sato, Asymptotically scale-invariant occupancy of 
phase space makes the entropy Sq extensive, cond-mat/0502274; T.S.Biro, G.Purcsel, 
G.Gyorgyi, A.Jakovac and Z.Schramd, Power-law tailed spectra from equilibrium; 
|nucl-th/0 510008, talk given at QM2005, Budapest (and references therein). 



5 




Figure 1: Examples of application of IT approach to different kind of data. Top-left: dN/dy 
for pp and pp collisions at different energies, fitted parameters are inelasticity and q = qi, 
cf., [HI for details. Top-middle: px distributions at different energies |7| (fitted parameter is 
q = q T which remains much smaller than q^. Top-right: dN/dy in the so called "tube model" 
of pA collisions, here K and R are fraction of energies of, respectively, incoming nucleon and 
"tube" used in hadronization process, cf., [S] for details. Bottom-left: p — (z/)nucleonic 
rapidity distributions for different values of hit nucleons in the tube model jS] . Notice that 
data on collisions of incoming proton with nuclear nucleon ( v = 1) disagree with data for 
pp collision at the same energy (lowest curve at upper-left panel) and cannot be therefore 
fitted - the probably reason is that "nucler nucleon" is a mixture of p and n. Bottom-middle: 
dN/dr] distributions in AA collisions [9 . Bottom-right: attempt to fit e + e~ data - as one 
can see they cannot be fitted properly by IT the method, cf. j^j. Here the best fit is for 
q < 1 which therefore limits available phase space. 
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